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ABSTRACT 

Assuming Coulomb-like as well as confining scalar potential, we have solved 
Shrodinger equation perturbatively in 1/toq with a heavy quark mass itiq. The 
lowest order equation is examined carefully. Mass levels are fitted with exper- 
imental data for D/B mesons at each level of perturbation. Meson wave func- 
tions obtained thereby can be used to calculate ordinary form factors as well as 
Isgur-Wise functions for semileptonic weak decays and other physical quantities. 
All the above calculations are expanded in I/tuq order by order to determine 
parameters as well as to compare with results of Heavy Quark Effective Theory. 

1. Introduction 

HQET (Heavy Quark Effective Theory), 1 is applied to many aspects of high en- 
ergy theories and many kinds of physical quantities of QCD which can be pertur- 
batively calculated in l/rriq. Especially those regarding to B meson physics, e.g., 
the Isgur-Wise function of semileptonic weak deacy processes B — > Dlv and the 
Kobayashi-Maskawa matrix element have been calculated by many people. How- 
ever, applications of HQET to higher order perturbative calculations are limited only 
to obtain forms of higher order operators, and their coefficients should be obtained so 
that results be fitted with experimental data. This is because most of the calculations 
based on HQET do not introduce heavy meson wave functions and hence there is no 
way to determine those coefficients within the model. 

In this paper, we would like to start from introducing phenomenological dynamics, 
i.e., assuming Coulomb-like vector and confining scalar potential to Qq bound states 
(heavy mesons), expand a hamiltonian in 1/mq then perturbatively solve Shrodinger 
equation in 1/rriQ. Angular part of the lowest order wave function is exactly solved. 
After extracting asymptotic forms of the lowest order wave function at both r — > 
and r — > oo and adopting the variational method, we numerically obtain radial part 
of the trial wave function which is expanded in polynomials of radial variable r. 
Then fitting the smallest eigenvalues of a hamiltonian with masses of D and D* 
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mesons, a strong coupling a s and other potential parameters are determined uniquely. 
Using parameters obtained this way, other mass levels are calculated and compared 
with experimental data for D/B mesons at each level of perturbation. Meson wave 
functions obtained thereby may be used to calculate ordinary form factors/Isgur-Wise 
functions for semileptonic weak decay processes and other physical quantities as well. 



2. Hamiltonian and Eigenvalue Equation 

Our hamiltonian is given by i 

H = (<V p q + f3 q m q ) + (a Q ■ p Q + f3 Q m Q ) + f3 q /3 Q S(r) 
+ j 1 - g [^q-^Q + (^q-ri)(dQ-n)]^V(r). 

When a heavy quark, Q, is treated as a non-relativistic particle, the hamiltonian is 
reduced into a 4 x 2 matrix operator and the Schrodinger equation becomes 

(#_ a + #„ + Hi + H 2 + ■ ■ •) 4> = M<j> = (m Q + E b )<p 
where m q — E b is a binding energy and scalar and vector potentials are given by 

S(r) = ^ + b, V(r) = ~^, and E b = E + ( E x + • - • . 
a z 3 r \m Q j 

and 

H-i = m Q 
H = a q -p + P q {m q + S) + V. 

Higher order terms like H\ and H 2 have some complicated stuructures and are not 
described here for simplicity. Subscripts i of Hi denote the order of 1/tuq. 

The -1st order eigenvalue equation is given by H^cpQ = mQ(f> , which is a trivial 
one. The next 0-th order becomes a non-trivial equation, H ^/^ m = E^ k - m , which 
can be reduced into 

fm q + S + V -d r + ^ \fu k (r)\ = fu k (r)\ 
V d r + t - mq -S + Vj\v k (r)J k \v k (r)J- 

This equation is numerically solved by taking into account the asymptotic behaviors 
at both r — > and r — ► oo and their forms are given by 

Uk(r),v k (r) ~ w(r) r 7 exp \ -{m q + b)r - ^ (J-^J j , 



where 7 = yk 2 — (4 a s /3) 2 and w{r) is some polynomial of r. In the above an 
eigenfunction, ^j m , can, in general, be given by 

•}„(?,«) =«.(r)rf„.(n) = i (_ iv ;;yi. n) ) K?„(n), 



where the function y^ m {Q) is a linear combination of Yj m and vector spherical har- 
monics a ■ Xjj(Q) and a ■ Yj? ± JQ), which satisfies 

— # 

with = ±j, or ± (j + 1) and £ is an angular momentum. The first few series of y| m 
functions are given by 

1 ,* * ..-i 1 



Vo o = nT (« • ^) > % o 



4^ v " yuu v^' 



where = -(ra • a) y) m , cr± = ia x ± ia y )/2, and n = f/r. 

Matrix elements of interaction terms among eigenfunctions are expressed as below. 
Degeneracy can be resolved by diagonal elements of Hi and H 2 with respect to the k 
quantum number. Inclusion of off-diagonal elements of H\ and H2 are absorbed into 
wave function corrections. Calculating all the matrix elements from the hamiltonian 
given above, total mass matrix is given by 

/ b_i+b_i, 00 \ 

e-i+E-i,! v-1,2 

Bi+Ei, 

e 1+ e 1a vi,_ 2 

V-2,1 B-2+S_ 2>1 



v 2 ,-i £ 2 +£2,i 

V 0000 00 e 2 +e 2 , 2 J 

where 

Zy 1 lyy, _J_ , __J_ T^ 2 ) F - _L_ T/ - T/W _L_ T/( 2 ) 

superscripts mean the order of 1/rriQ, k and k! stand for k quantum number, and j 
for the total angular momentum. 

Corrections to wavefunctions and eigenvalues can be calculated by applying the 
ordinary perturbation. For instance, in the above 8x8 mass matrix up to 0(l/m c ), 
the wave function for (2, 2) element is given by 

* 2 = M>r>, ^) - 7.92 x 10-^ 2 lm (r,n) 

m,, 



= *r;(^)-i.i9xp*L(^), 



where 



a 2 

^- = 0.299, a = 2.10 GeV" 1 , b = -0.084 GeV, 
4% 



m u = 0.01 GeV, m c = 1.461 GeV, m b = 4.92 GeV. 

If D and D* are identified as ^/qo 1 an d ^im respectively, D and B meson masses are 
calculated to be 

m D = 1.867 GeV, m D * = 1.9241 GeV, m B = 5.279 GeV, m B * = 5.298 GeV. 

3. Comments and Discussions 

One can easily see degeneracy among the lowest lying pseudoscalar and vector 
states as follows. Define 

\P) = |D± D°) = %l \V, m) = \D*) = Vjl, 

where is an eigenfunction obtained in the last chapter and quantum number k 
can take only ±j, or ± (j + 1). Since these states have the same quantum number 
k — — 1, these have the same masses as well as the same wave functions up to the 0-th 
order calculation in I/tuq. One needs to develop perturbation of energy and wave 
function for each state in terms of 1/rriQ to obtain higher order corrections. 

Next we would like to discuss qualitative features of form factors functions. Let 
us think about to calculate form factors for semileptonic decay of B — > Dtv. Taking 
a simple form for the lowest lying wave function both for B and D as 

q?™ ~ e~ b2 r2/2 , 

where a parameter b is determined by a variational principle, 5(^ 1S ^H^/ 13 ) = 0. 
Then form factors/Isgure-Wise functions are given by 



F(q ) ~ exp const. m q (q - g max ) , or £(u) ~ exp 



const. m 2 q (1 — uj) 



where 

<? 2 = (Pb ~ Pd) 2 , uj = v B ■ v D , g^ ax = (m B - m D f <-> ^ max = 1. 

This means behavior of form factors strongly depends on a value of light quark mass 
m q . One needs current quark mass for m q to reproduce heavy meson mass spectrum, 
while constituent quark mass is apparently used to calculate form factors in all pub- 
lished papers. In order to explain this situation, one is forced to use running mass i 
m q (q 2 ) for light quark mass. 
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